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1 Introduction 

Quantum Chroniodynamics (QCD) is nowadays settled as the theory of strong interac- 
tions within the Standard Model (SM) of the elementary particles. To a large extent, this 
achievement is the result of the theoretical and experimental progress in hadronic physics at 
high energies. In this kinematic regime, QCD is synonymous of 'perturbative QCD'. This 
talk will cover a selection of topics, rather than results, in perturbative QCD. More com- 
prehensive reviews of recent results on jet physics □, unpolarizedH and polarizedu nucleon 
structure functions, photon structurea, diffractive processesH and heavy-quarkQ production 
and decay are presented in other contributions to this Symposium. 

Calculations at the lowest order (LO) in QCD perturbation theory give only the order 
of magnitude of hard-scattering cross-sections. The theoretical accuracy of perturbative- 
QCD predictions is instead controlled by the size of the next-to-leading order (NLO) and, 
in general, higher-order contributions. Section |2| is devoted to a review of the theoretical 
progress in perturbative calculations and also includes the results of recent measurements 
of as- 

Soft-gluon radiation is a source of large higher-order corrections for hard-scattering 
processes near the exclusive boundary of the phase space. In these cases, summation of 
the corrections to all orders can be important to improve the accuracy of the perturbative 
approach. In Sect. ^, some predictions based on soft-gluon resummation are briefly reviewed 
and their impact on the analysis of experimental data from high-energy colliders is discussed. 

Low-x physics is a topic at the border of hard and soft hadronic interactions. Section ^ 
summarizes present analyses of the low- a; behaviour of the nucleon structure function and 
outlines recent theoretical developments in the BFKL formulation of small- a; dynamics. 

At high energies, non-perturbative phenomena affect perturbative predictions by contri- 
butions that are suppressed by inverse powers of the energy. Some recent theoretical ideas 
and phenomenological studies to quantify power corrections are discussed in Sect. ^. 

A world summary of as determinations is presented in Sect. ^ and some concluding 
remarks are left to Sect. 0. 



2 Higher-order calculations in perturbation theory 

The evaluation of perturbative corrections to hard-scattering cross-sections requires the 
computation of higher-order Feynnian diagrams that involve real and virtual partons. In this 
computation one has to deal with different kinds of singularities. The customary ultraviolet 
singularities, present in the virtual contributions, are removed by renormalization. By 
adding real and virtual terms, the infrared divergences cancel in inclusive cross sections, 
while the left-over collinear singularities are factorized in the process-independent parton 
distributions and fragmentation functions. 

Because of this complicated pattern of singularities, it is natural to divide QCD observ- 
ables into two different classes, according to their degree of inclusiveness. 

2.1 Completely inclusive quantities: NNLO predictions and recent as determinations 

Fully inclusive quantities are infrared- and collinear-safe observables that depend on a single 
momentum scale. They can be expressed as a simple power series expansion in as- 

The best known observables of this type are the total hadrouic cross-section in e'^e~ 
annihilation and the hadronic branching ratio of the Z'^ . In Ref. Q, these observables were 
computed up to next-to-next-to-leading order (NNLO) in perturbation theory, i.e. to rel- 
ative accuracy 0{a^g) with respect to the lowest-order approximation. These calculations 
have recently been confirmed by an independent re-evaluation B. Using these results, one 
can parametrize B the hadronic branching ratio of the Z^ as follows 
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where the factor Rq includes the electroweak radiative corrections and, in particular, de- 
pends on the values of the masses of the top quark and Higgs boson. The term 0{{\/MzY) 
stands for non-perturbative power corrections (cf. Sect. ^). 

Other quantities, which have been computed at NNLO, are the hadronic width of the r 
leptonli3'EJI and the deep inelastic scattering (DIS) sum rulesEJ, namely the Gross-Llewellyn 
Smith sum rulellj and the polarized llil and unpolarized Bjorken sum rules. 

Fully inclusive observables are the simplest quantities that can accurately be evaluated 
in QCD perturbation theory. Their (relative) simplicity from the computational viewpoint 
follows from kinematics. Since these observables are completely inclusive, no phase-space 
restriction has to be applied. Real and virtual contributions can be combined at the inte- 
grand level and this produces the cancellation of infrared and collinear singularities before 
performing the relevant phase-space integrations. Owing to these features, general tech- 
niques have been available for some time Ej to carry out NNLO calculations in anaWtic 
form. These techniques are suitable for automatic implementations in computer codes H^. 

The high accuracy of the theoretical predictions for these observables derives from the 
availability of NNLO calculations and from the fact that non-perturbative power corrections 
can be controlled by operator product expansions 113 (OPEs) and are relatively small (in 
general, p>2 and, in particular, p = 4 in Eq. @)). Because of these reasons, fully inclusive 
observables are particularly suitable for as determination. 



The cleanest and, in principle, most accurate determination of as is that performed at 
the Z^ peak. Here, effects of power corrections are strongly suppressed because of the large 
value of the Z mass, and the theoretical uncertainty due to perturbative-QCD contributions 
beyond NNLO is estimated to be {Aas)QCD = ±0.002. In practice, however, the value of 
as obtained in this way is quite sensitive to the assumption that electroweak interactions 
are accurately described by the SM. For instance, extracting as{Mz) from Eq. (||), one 
hasta 

as = {as)s.M. - 3.2 — ^ . (2) 

[I- ha.d)S.M. 

One can see that a variation of the hadronic width by one per mille with respect to the value 
expected from the SM produces an effect on as that is larger than the QCD uncertainty. 

Using a fixed top mass, mt = 175.6ib5.5GeV, and the experimental value of Rz reported 
at the 1997 summer conferences 113, one obtains from Eq. (|l]) 

asiMz) = 0.124 ± 0.004 (exp.) ± 0.002 (m^) , (3) 

where the central value corresponds to the Higgs mass rriH = 300 GeV and the second error 
is due to the variation of mn in the range 60 GeV < ttih < 1 TeV. The sensitivity of as 
to the SM assumption can be reduced by considering a global fit to electroweak data. In 
this case as mainly depends on Rz, the total width Tz and the peak value of the hadronic 
cross section a^. The simultaneous fit of 'mt,as,mH givesEJ nit = 173.1 ± 5.4 GeV and 

asiMz) = 0.120 ± 0.003 (exp.) , mn = 1151^66 GeV . (4) 

The difference between the values (^) and (^) for as comes from a shift (~ —0.002) produced 
by the new entry Tz and a further shift (~ —0.002) due to the different central value of nin- 
This shows that, within a global SM fit, as is still quite sensitive to uih- This dependence 
can be parametrized as follows Ii3 

{Aas)niH = 0-0023 xh(1 + 0.2 xh) , xh = ln(mH/100GeV) . (5) 

Much theoretical workc3 has recently been devoted to the calculation of quark-mass cor- 
rections to current-current correlators. These corrections are important for measurements 
of as from low-energy e^e~ data. 

These results and estimates of non-perturbaiive contributions f along the lines of similar 
analyses for the hadronic decay of the r leptont^) have been usedcil to determine as from 
the e~^e~ hadronic cross-section at the T resonance. This NLO determination gives 

a5(4.1GeV) = 0.2281°;°^^ , (6) 

where the error is dominated by the theoretical uncertainty. The corresponding value 
evolved to the Z'^jnass is as{Mz) = 0.119lo!oo8 ^^*^ ^^^ ^^ ^^ compared with a previ- 
ous determination E3, as{Mz) = 0.109 ± 0.001, based on LO predictions. 

Measurements of as from the continuum in e^e~ annihilation at low energies have 
typically large experimental errors due to poor statistics. The new measurement submitted 
by the CLEO Collaboration to this Symposium isc3 

05(10.52 GeV) = 0.20 ± 0.01(stat.) ± 0.06(syst.) . (7) 



Note that, because of the large amount of data collected at CESR just below the T res- 
onance, the statistical error is no longer dominant. The systematic error is theoretical 
(uncertainties in QED radiative corrections) and experimental (estimates of backgrounds 
and detector efficiencies). 

2.2 Inclusive quantities: NLO calculations and general algorithms 

QCD calculations beyond LO for inclusive quantities are much more involved. Owing to the 
complicated phase space for multiparton configurations, analytic calculations are in practice 
impossible for most of the distributions. Moreover, infrared and collinear singularities, 
separately present in the real and virtual contributions at the intermediate steps, have to 
be ffist regularized by analytic continuation in a number of space-time dimensions d = 4 — 2e 
different from four. This analytic continuation prevents a straightforward implementation 
of numerical integration techniques. 

Methods to overcome these problems are known. They consist in using hybrid analyt- 
ical/numerical procedures: one must somehow extract and simplify the singular parts and 
treat them analytically; the remainder is treated numerically, independently of the full com- 
plications of the actual calculation. These methods were first ^3 used to evaluate 3-jet cross 
sections in e"*"e~ annihilation and were then applied to other cross sections c3, adapting the 
method each time to the particular process. This is very time-consuming and has required 
lot of effort to produce new NLO calculations. 

Only recently has it become clear that these methods are generalizable in a process- 
independent manner. The key observation is that the singular parts of the QCD matrix 
elements can be singled out in a general way by using the factorization properties of soft 
and collinear radiation. Owing to this universality, the methods have led to general algo- 
rithms c3'ElllHa for NLO QCD calculations. 

The various algorithms use different methods and techniques (phase-space slicing E3, 
subtraction method E3, dipole formalism IH2) to achieve a common goal. All the analytical 
work that is necessary to evaluate and cancel the infrared singularities is carried out once 
and for all. The final output of the algorithms is given in terms of effective matrix elements 
that can be automatically constructed starting from the original (process-dependent) ma- 
trix elements and universal (process-independent) factors. The effective matrix elements 
can be integrated numerically or analytically (whenever possible) over the available phase 
space in four dimensions to compute the actual value of the NLO cross section. If the nu- 
merical approach is chosen, Monte Carlo integration techniques can be easily implemented 
to provide general-purpose Monte Carlo programs for carrying out NLO calculations in any 
given process. 

Using these algorithms, the computation of inclusive quantities at NLO essentiallY 
amounts to the evaluation of the original matrix elements. Since efficient techniques c3 
(based on helicity formalism and colour-subamplitude decompositions) are available for 
calculating real matrix elements, the computation of the matrix elements for the virtual 
contribution remains the only real obstacle to perform new NLO calculations. 

In recent years, this obstacle was greatly reduced by the introduction of new tools c3, 
inspired by string-theory methods, for the evaluation of one-loop amplitudes. One-loop 



matrix elements involving up to five massless partons are known O and the complication of 
those with four massless partons and a vector boson has been completed recently c3. 

The relevance of the theoretical progress that I have briefly outlined is witnessed by 
the accelerated production rate of new NLO calculations. These include, for instance, 4- 
jet cross sectionsE3 and mass quark corrections^ to 3-jet observables in e"'"e~ annihilation. 
The complete calculation for 3-jet cross sections in hadron collisionsE^ is expected to appear 
soon. 

The theoretical accuracy of NLO predictions for inclusive observables can reach the 
level of 10-20%, although in some notable cases (e.g. production cross sections of charm 
and bottom quarksQ and direct photonsEl) theoretical uncertainties are much larger and the 
agreement with data is poor. This accuracy is very important for physics studies within and 
beyond the SM. The comparison between NLO predictions and high-energy collider data is 
discussed in Ref.u and in the following sections. The new calculations mentioned above are 
extremely valuable for further studies such as, for instance, to obtain stronger constraints 
on the elusive light gluinoE3c3 and to measure the running of the 6-quark massoO. 

The extension to higher orders of the NLO techniques described in this section is still 
a challenge for theorists. The high experimental accuracy of LEP and SLC data demands 
improved perturbative-QCD predictions and efforts in this direction. 

3 All-order resummation, high Et, high x 

Higher-order perturbative contributions are certainly important for studies of QCD observ- 
ables close to the exclusive boundary of the phase space. Although infrared and collinear 
singularities cancel in inclusive cross sections upon adding real and virtual contributions, 
in this kinematic regime real emission is strongly inhibited. The ensuing mismatch of 
real and virtual corrections generates double-logarithmic terms of the type {asL'^)"', where 
L = In 1/y, and y generically denotes the distance from the exclusive boundary. For in- 
stance, y = ycut can be the resolution parameter for the transverse size of jets in e~^e~ 
annihilation, or y = 1 — 2Et/vS, where Et is the transverse energy of jets produced in 
hadron collisions at the centre-of-mass energy v5, ov y = 1 — x, where x is the Bjorken 
variable in DIS processes. 

In all these processes, when y <^ 1 the logarithmically-enhanced terms of infrared origin 
spoil the convergence of the fixed-order expansion in as- Accurate predictions require the 
evaluation and (whenever feasible) the resummation of this class of contributions to all 
orders in perturbation theory E3. In the following I briefly discuss some examples of soft- 
gluon resummation in lepton and hadron collisions. Transverse-momentum distributions of 
vector bosons and di-photon systems are reviewed elsewhere EJ in these proceedings. 

3.1 Jet rates and event shapes in e+e^ annihilation 

A detailed understanding of logarithmically-enhanced terms exists for jet rates E3 and for 
many shape variables e3'e2I in e^e~ annihilation in the two-jet limit. In these cases all-order 
resummation takes an exponentiated form. For instance, using the /cx-algorithm Ej'e3 to 



define jets, the 2-jet rate R2 is given byEja 



R2[as{Q),ycut) = ^^ = e^P< / -- as{k)ln—^ + ...} , (8) 

o"tot I TT Jq^/y^ k k^ I 



where ycut is the jet resolution parameter, Q is the e^e~ centre-of-mass energy and the dots 
denote subleading contributions. Note that the exponent involves a double-logarithmically- 
weighted integral of the QCD coupling. The effective scale k at which as is evaluated corre- 
sponds to the transverse momentum exchanged at the elementary QCD vertices rather than 
to the total momentum Q. Using the QCD running, as{k) = as{Q)/{'i- + 2PQas{Q) Ink/Q), 
one can nonetheless perform the integral in Eq. (|8|) and explicitly obtain R2 as a function 
of as{Q) and L = In 1/ycut- 

In general, shape variable and jet rates can written as follows o 

R{as, y) = C{as)nas. L) + D{as. y) , (9) 

00 
^{as,y) = e^Y>{Lgi{asL)+g2{asL) + ---] , C(a5) = 1 + ^ C„ag , (10) 

ra=l 

where D{as,y) vanishes as y — > order-by-order in perturbation theory. In the expression 
@ the singular In y-dependence is entirely included in the effective form factor S, given in 
Eq. (p!o[). The function L gi resums all the leading contributions a'gL^^^, while g2 contains 
the next-to-leading logarithmic terms a^L^ and so forth. Equation (^) represents an 
improved perturbative expansion in the two-jet region y ^ 1. Once the functions gi have 
been computed, one has a systematic perturbative treatment of the shape distribution 
throughout the region of y in which asL ^ 1, which is much larger than the domain a^L^ <C 
1 in which the fixed-order expansion in as is applicable. Furthermore, the resummed 



expressions (^) and (10) can consistently be matched with fixed-order calculations. In 
particular, one can consider the next-to- leading logarithmic approximation (NLLA) as given 
by the functions gi and g2 and combine them with the complete NLO results r tt^I (after 
subtracting the resummed logarithmic terms in order to avoid double counting), to obtain 
a prediction (C'(q;|)-|-NLLA) which is everywhere at least as good as the fixed-order result, 
and much better as y becomes srnall. 

Detailed experimental studiese^ performed at LEP 1 and the SLC proved that resummed 
predictions substantially improve the QCD description of the data. These analyses led 
to determinations of as from hadronic final states with reduced theoretical uncertainty 
with respect to pure NLO calculations. The combined value from SLD and the four LEP 
experiments is as{Mz) = 0.122 ± 0.001(exp.) ± 0.006(th.). 

Resummed calculations were also compared with e^e~ data at lower energies, namely 
at PEP§ (Q = 29 GeV) and TRISTAN§ (Q = 58 GeV). A re-analysis of the data recorded 
in 1981-1986 by the JADE detector at PETRA has been presented recently eHI. This anal- 
ysis provides measurements of as at three different centre-of-mass energies, 05(22 GeV) = 
0.161, a5'(35GeV) = 0.143, 0^(44 GeV) = 0.137. The energy dependence of these values is 
in agreement with the QCD expectation and corresponds to 

as{Mz) = 0.122t0:00^ , (11) 
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Figure 1: The differential 2-jet rates, di?2/rfy23 (1/23 = J/cut), measured at Q = 44,35 GeV, are shown after 
correction to the parton level. The solid and dashed lines correspond to the results of fits with resummed 
and fixed-order calculations, respectively. Note the extension of the fit range (denoted by the arrow) towards 
the sniall-y23 region in the case of resummed predictions. 





Thrust (T) Heavy Jet Mass (p) 

Figure 2: Measured distributions of thrust T (left) and scaled heavy jet mass p (right) in comparison with 
C(a|)+NLLA QCD predictions at Q = 172 GeV. 



where the error is dominated by the theoretical uncertainty. The result ( |llD improves a 
previous determination, as{Mz) = 0.119 it 0.014, based on NLO predictions. The improve- 
ment is due to the use of more observables (see, e.g., Fig. |l|) and to the inclusion of more 
detailed systematic studies. 

The most recent data E^ on hadronic events from LEP 2 around Q = 133, 161 and 
172 GeV are well described by QCD predictions. Resummed calculations are quite valuable 
because they can be used for studies of hadronic observables close to the 2-jet region (Fig. |2|), 
thus increasing the (otherwise very limited) statistics. The ensuing first measurements of 
as at these new energies are summarized in Sect. |^. Assuming the QCD running, the 
combined LEP 2 average gives asiMz) = 0.115 ± 0.008. 



3.2 Hadron collisions at high Et,Q^ and x 

At high-energy hadron colliders, the region close 
space is the most sensitive to possible signals E 
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Figure 3; One-jet inclusive cross-section in pp collisions at 1.8 TeV. 



dictions are thus important not only to test the SM but also to correctly estimate SM 
backgrounds. Note that these predictions depend on an extra non-perturbative input with 
respect to e^e~ collisions: the parton densities of the incoming hadrons to be convoluted 
with hard-scattering cross-sections. The present knowledge of the parton densities is re- 
viewed elsewhere □ in these proceedings. 

The top quark is produced relatively close to threshold at the TevatronH. Calculations 
of its production cross-section, based on NLO QCDEil and including soft-gluon resummation, 
were performed in Refs.EJ~Kj. The results disagree on soft-gluon effects. The disagreement 
does not regard the general framework e30 to carry out the resummation, but the proper 
way to implement resummed formulae in actual computations of hadronic cross-sections. 

The differences between Refs. E3'E3'E3 do not have a large impact on top-quark phe- 
nomenology: the various numerical results are consistent within uncertainty estimates and 
no sizeable reduction of the present experimental error is expected. However, the use of 
different approaches EJ'E3'l3 to soft-gluon resummation in hadron collisions can be more rel- 
evant in other processes as, for instance, in the case of jet production at large transverse 
energy Et- 

The one-jet inclusive cross-section measured at the TevatronEI impressively agrees with 
NLO QCD predictions!^ over almost 9 orders of magnitude (Fig. ^). However, the CDF 
Collaboration reported u an excess of events at high Et (-Et ^ 250GeV). The present 
situationlllEj can be summarized as follows (Fig. ^): CDF still sees an excess while DO does 
not and, nonetheless, CDF and DO data are consistent within experimental errors (Fig. |5|). 
Previous data from the two experiments could not be compared directly as they probe 
different regions of jet pseudorapidity tj. A new DO analysis in the same pseudorapidity 
range as that of CDF shows (Fig. |5|) that the CDF data lie above the DO results, but within 
the DO uncertainty band. 

The origin of the differences between these experimental results has to be further clar- 
ified. The angular distribution for two-jet events Ej is in good agreement with NLO QCD 
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Figure 4: Relative difference between the inclusive jet cross-section and NLO predictions. CDF (left) and DO 
(right) data refer to the pseudorapidity ranges 0.1 < \ri\ < 0.7 and |7;| < 0.5, respectively. CDF systematic 
errors are not shown. 

calculations, suggesting that the possible high-E'y excess may have an explanation within 
the framework of QCD rather than originating from new physics. Independently, the accu- 
racy of the QCD predictions has to be carefully estimated. 

Using different conventional sets E3'E3cZI of parton distributions and different values of 
factorization/renormalization scales, NLO predictions for the one-jet cross-section in the 
range 50 ^ Et ^ 400 vary by approximately 15%. As for sofLgluon effects, one can at 
present estimate an increase of the cross section by ~ 20% e3 or by less than 10% e3, 
depending on the resummation approch used. Further investigations are warranted. The 
fact that the estimate in Ref.l^ is within the variation of the NLO predictions may suggest 
that higher-order QCD corrections are under control. 

Note, however, that the true uncertainty due to parton densities can be larger. As shown 
by the CTEQ Collaboration^, including the CDF data in global fits to parton distributions, 
there is enough flexibility to increase the gluon density at large x and enhance by 25-30% 
the NLO predictions for the single-jet distribution at high Et- 

The HERA experiments have reported e3 an excess of DIS events at large values of 
Q2 (15 000^Q2^50 000GeV2) and x (0.5^x^0.7). This has raised considerable inter- 
est in the high-energy physics community!^. The measured excess has decreased as the 
integrated luminosity increased Lj, but it still deserves attention Ea. 

The observed excess is with respect to the SM expectation based on QCD predictions. 
Their accuracy relies on estimates of beyond-NLO perturbative corrections and on accurate 
determinations of non-perturbative parton densities. 

In DIS processes at large x, logarithmically-enhanced contributions due to soft-gluon 
emission can be sizeable. We have implemented Ell the corresponding resummation formu- 
lae E3 to estimate the order of magnitude of the effect in the relevant HERA region. The 
results of this excercise are illustrated in Fig. ^. The overall sign of the resummation effect 



9 



0.2 - 



„v 



nna Oi 



Band shows D0 systematic error ~~~~ 

-OA - 

_ra 100 150 200 2510 300 350 400 450 500 

E^(GeV) 

Figure 5; Relative difference between CDF data and a fit of the DO data in the range 0.1 < \r]\ < 0.7. The 
systematic band is mainly due to the DO jet energy scale uncertainty. 

on the measured structure functions is negative because multigluon radiation increases the 
scaUng violation at large x. The decrease in the predictions is nonetheless extremely small 
for a;;^0.7, so that the NLO evolution is safe in this kinematics region. 

The main source of uncertainty in the SM predictions^ thus comea^from the extraction 
of the parton (quark) densities from low-Q^ data. From a dedicated fitE3 to SLAC, BCDMS, 
NMC data at large x, the uncertainty is estimated to be ~ 10%. This value is larger than 
the differences (~ 6%) obtained by calculations that use different setsE3£3E3 of conventional 
parton densitiesfcf. Fig. 3 in Ref. O). Sizeable additional effects due, for instance, to^ 
new componentll^ in the quark densities at extremely large values of x, can be excluded E3 
by combining HERA data in the neutral- and charged-current channels. 

4 Low-x physics and structure functions 

Strong-interaction dynamics in the small-x regime l3 is a main challenge to QCD. Much 
progress in the field has been prompted by measurements of stucture functionsH and diffrac- 
tive interactions u at HERA in a kinematic range that extends down to x ~ 10"'^. The 
present status of small-x physics for hard-scattering processes involving two transverse- 
momentum scales (e.g. forward jet production in DIS, two-jet inclusive mesurements at 
large rapidity gaps) is reviewed in Ref. EJ. 

I shall consider single-scale processes and, in particular, DIS structure functions at low 
X. In this case, one can studyH the transition between the perturbative and non-perturbative 
regimes, and this may eventually lead to a QCD understanding of soft hadronic physics at 
high energy. 

4.1 Confronting DGLAP with BFKL 

The QCD analysis of low-x structure functions has to deal with a non-trivial interplay 
between perturbative evolution towards the hard scale Q^ and non-perturbative behaviour 

10 



Figure 6: The ratio F'"'''{x,Q'^)/F^^°{x,Q'^) at Q^ = 30 000 GeV^ The DIS structure functions F^^'^ 
and F'"^'^ are respectively obtained by Q^-evolution in NLO and including soft-gluon resummation. The 
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of the parton densities. This comphcates the answer to a basic question that regards the 



region of applicability of the QCD parton model. How far in x can we safely use 



evolution equations 



he DGLAP 

effects? 



;II3 before they have to be supplemented with BFKL-type 
To discuss this point let me recall the main steps in the QCD analysis of the structure 
functions. The measurement of the proton structure function ^2(2;, Q"^) ~ fs{x^ Q^) directly 
determines the sea-quark density fs = x{q + q). Then one can use the DGLAP evolution 
equations (the symbol denotes the convolution integral with respect to x), 



dfgix, Q^)/dln Q2 ^ Pg^ 55 /^ + Pg^ 55 fg , 



(12) 
(13) 



to extract a gluon density fg{x, Q^) = xg{x, Q^) that agrees with the measured scaling vio- 
lation in dF2{x, Q'^)/d\\iQ^ (according to Eq. (|l^)) and fulfils the self-consistency equation 

(HD- 

The perturbative-QCD ingredients in this analysis are the Altarelli-Parisi splitting func- 
tions Pab{oi.s{Q)-,x). They are computable as power-series expansions in as and are known 
up to NLO accuracy. The truncation of the sphtting functions at a fixed perturbative or- 
der is equivalent to assuming that the dominant dynamical mechanism leading to scaling 
violations is the evolution of parton cascades with strongly-ordered transverse momenta. 
However, at high energy this evolution takes place over large rapidity intervals (Ay ~ In 1/x) 
and diffusion in transverse momentum becomes relevant. Formally, this implies that higher- 
order corrections to Pab(as,x) are logarithmically enhanced: 



Pab{oi.s,x] 



^ ^ + ^ (05 In x) + . . . + ^ (05 In x)" + . . . 

XX X 



(14) 



At asymptotically small values of x, resummation of these corrections is mandatory to 
obtain reliable predictions. 

Small-x resummation is, in general, accomplished by the BFKL equation O, whose 
structure is completely known only to leading logarithmic (LL) accuracy. In the context of 
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structure-function calculations, the BFKL equation provides us with resummed expressions 
for the splitting functions. 

In the small-x region the gluon channel dominates. Considering the fixed-order expan- 
sion of the splitting function, one-gluon exchange gives Pgg{as,x) ~ as/x. Then, assuming 



a flat x-behaviour of fg at low momentum scales, the evolution equation (|13|) produces_a 
gluon density fg ~ ex.p{^/\n 1/x) that steeply increases in 1/x at higher momentum scalesE^. 
The steep behaviour drives strong scaling violations in F2 and leads to structure functions 
that strongly rise as x decreases. This prediction is usually referred to as DGLAP prediction. 

The increase of the gluon density is steeper after BFKL resummation. Summing the LL 
terms in Pgg, one eventually gets the power-like behaviour fg ~ x"'^", where A^ ~ 2.65 a^. 

The fixed-order approach has been extensively compared with structure function data 
during the last few years. The NLO approximation of the DGLAP equations describes very 
wenHililll-iil the HERA data (Fig. 0), down to low values of Q^ ^ 2GeV2. The NLO 
QCD fits simply require slightly steep input parton densities at these low momentum scales. 

From a phenomenological viewpoint, one may conclude that there is no need for BFKL- 
type corrections to scaling violations in the HERA kinematic region. The main reason for 
this is that the LL power-like behaviour of the BFKL gluon density is valid in the asymptotic 
regime and the approach to asymptotia is much delayedE3~H3 by cancellations of logarithmic 
corrections that occur at the first perturbative orders in the gluon splitting function Pgg. 

Nonetheless, a better way to estimate the relevance of BFKL-type corrections is to 
quantify the theoretical uncertainty of the NLO predictions, for instance by comparing LO 
and NLO results. Doing that, one can thus argue that this uncertainty is sizeable. This 
feature has recently begn re-emphasized in Ref.Ell, but it was evident since the still-successful 
GRV parametrizationS of parton densities. Going from LO to NLO, one can obtain stable 
predictions for the proton structure function F2, but one has to vary the parton densities a 
lot, in particular the gluon. As shown in Fig. ^, the NLO gluon density sizeably differs from 
its LO parametrization, not only in absolute normalization but also in x-shape. This can be 
understoodl^B from the fact that the scaling violation of F2 is produced by the convolution 
Pgg fg (sBB thc rlghf-haud side of Eq. (0)). The quark splitting function Pgg behaves as 



Pqg{as,x):^asP^^\x) 



1 + 2.2^^1 + ... 

TT X 



(15) 



where the LO term Pgg (x) is flat at small x, whereas the NLO correction is steep. To 
obtain a stable evolution of F2, the NLO steepness of Pgg has to be compensated by a gluon 
density that is less steep at NLO than at LO. 

In the large-x region, there is a well known correlation between as and fg. At small 
X, there is an analogous strong correlation between the x-shapes of Pgg and fg. In the 
fixed-order DGLAP analysis of F2, large NLO perturbative corrections at small x can be 
balanced by the extreme flexibility of parton density parametrizations. 

This has to be kept in mind when concluding on the importance of BFKL dynamics. The 
NLO steepness of Pgg is the lowest-order manifestation of next-to-leading BFKL corrections 
in the quark channel. Using A:_|_-facto 



to all orders E2l and implemented 



rizationEll methods, these corrections were resummed 
|89|-p2| jj^ studies of structure functions. The results of 
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Figure 7: NLO QCD fit to low- a; structure function data for difTerent values of Q 
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Figure 8: Comparison between the LO (GRV94LO) and NLO (GRV94HO) GRV parametrizations of the 
gluon and sea-quark densities at Q^ = 5 GeV^ . 
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Figure 9; Q^ dependence of the proton structure function F2{x,Q^) at various values of x. Data from 
the E665, HI and ZEUS experiments are compared with a QCD fit that includes resummation of small- a; 
contributions to LL accuracy in the gluon channel and to NLL accuracy in the quark channel. 
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a recent fit E3 are shown in Fig. 0. The present overall picture is that small-x resummation 
leads to a description of the HERA data that is comparable to that from NLO analyses, 
provided that the input gluon density is less steep than the NLO fg. 

More accurate data on F2 are necessary to clarify the phenomenological relevance of 
BFKL effects. Measuremeiits of other inclusive observables (for instance, the longitu- 
dinal structure function Fl Ej'o) can be valuable to disentangle perturbative from non- 
perturbative dynamics, that is, to reduce the correlation between Pgg and fg in F2 analyses. 
At the same time, theoretical progress in higher-order calculations is very important. 

4-2 BFKL dynamics at NLL accuracy: recent theoretical progress 

The complete evaluation of BFKL contributions to next-to-leading logarithmic (NLL) accu- 
racy requires the computation of basic building blocks from tree-level, one-loop and two-loop 
amplitudes at high energy. This difficult calculational program l3 t3, led by the work of 
Fadin and Lipatov, has been completed!^. 

We are now in a position to start detailed investigations on phenomenological and 
conceptual issues related to the resummation of high-energy logarithms. 

The computation of next-to-leading correction .s should permit a quantitative control 



of normalizations, scales and factorization-scheme E3't23't22l uncertainties. A preliminary 
evaluation of the NLL corrections to the power-like behaviour of the gluon density has been 
performed in Ref.". The effective power Xg of the asymptotic x-shape of fg can be written 
as Xg = 2.65 as{Q){l — c as{Q)) and, considering the scale-invariant part of the NLL BFKL 
contributions, the correction to Xg is negative and estimated to be quitelaxgeH (c ~ 3.5). 



Negative corrections to Xg are also obtained by related investigations E23't23 on coherence 



effects liiia and the CCFM equation 
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A general theoretical understanding of NLL terms can clarify to what extent the BFKL 
formalism can be used in a purely perturbative framework and how it can be matched 
to the non-perturbative infrared regime. In the context of structure-function studies, one 
should thus be able to precisely identify the kinematic range in x and Q^ where pRrt urhatj ve 
factorization is valid and not spoiled by A;_L-diffusion in the low-momentum region 
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5 Power corrections 

A generic infrared- and collinear-safe observable R that depends on some large momentum 
scale Q has the following expression 

RiQ) = RpertiasiQ)) + Rnorr-pcrtiQ) ■ (16) 

The term Rpcrt denotes the perturbative component that can be calculated as power-series 
expansion in as{Q) and, thus, behaves as (l/lnQ)"- The remaining contribution i?non-pert 
is due to non-perturbative phenomena (hadronization, multiparton scattering, ...) and is 
power-behaved, i.e. proportional to {1/Qy. Since the power p is positive, i?non-pert is 
suppressed when Q — > oo but it can be quantitatively relevant at finite values of Q. 
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We have precise theoretical information on^jion-pei-t only for the processes in which 
OPEs are validO. In these cases, one can writell23 

Rnon-peAQ) ^ Y. ^p ( ^Y ^^P^ ' ^^^^ 

p>2 ^"^^ 

where (Op) denotes the vacuum expectation value of some local operator Op, obtained by 
the basic quark and gluon fields. Note that in Eq. ( p!7| ) we have p > 2, because in the theory 
there is no (quasi-)local operator with dimension smaller than 2. Note also that OPEs apply 
to few quantities such as the total cross section in e^e~ annihilation (p > 4) and the DIS 
structure functions {p > 2). For all the other quantities, power-correction contributions 
are usually estimated by using phenomenological approaches (e.g. hadronization models of 
Monte Carlo event generators in e~^e~ annihilation). 

It is evident that we need a better general understanding of power-suppressed contri- 
butions. On theoretical grounds, one would like to justify power corrections with p < 2, 
which are experimentally observed in e^e^ data on event shapes (see Fig. p^ ), and possibly 
evaluate their size. 

A handle on i?non-pert can be provided by the study of i?pcrt- It is known that the 
a5-series for the pertuibative component i?pert is not a convergent series but rather an 
asymptotic expansion t23. This implies that it is not defined in an unambiguous way. Since 
R is a physical observable, any ambiguity in i?pert has to be cancelled by a corresponding 
ambiguity in i?non-pert- Examining the series for -Rpert at large perturbative orders, one can 
thus extract some information on i?non-port- 

Two known sources of non-convergent behaviour of the perturbative expansion are in- 
stantons and (infrared) renormalons. Roughly speaking, 'instantons' are related to large 
perturbative coefficients due to the large number of Feynman diagrams at high orders. They 
lead to non-perturbative corrections that are strongly power-suppressed, i.e. that have high 
values of p. Infrared renormalons are related to the low- momentum behaviour of the running 
coupling as in Feynman diagrams. 

5.1 Infrared renormalons 



lot 



Much theoretical activity has recently been devoted to studying infrared renormalons 
A comprehensive review of the field is not feasible because even a simple list of the relevant 
references cannot fit into several pages of these proceedings. I shall limit myself to recalling 
few general points and phenomenological results. 

The basic idea of the renormalon approach to power corrections is the following. The 
evaluation of the lowest-order perturbative contribution to R{Q) amounts to integrating 
tree-level Feynman diagrams over the momentum k exchanged at the elementary QCD 
vertices (see the left-hand side of Eq. (|l^)). Higher-order contributions to R{Q) are given 
by more complicated Feynman diagrams, including those that produce the running of the 
coupling constant as- Therefore, one can approximate the effects of higher orders by the 
replacement as — > as{k) in the lowest-order term, as in Eq. (|l^): 

dk/ky ^ rQdk/ky 



RpertiasiQ)) r^ as 1^ — (t^) +••• high^rders /„ 77 ( 7T ) ^s(.k) . (18) 
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The perturbative series Rpert{cts) = Sn -^nog is then obtained by inserting the QCD ex- 
pression for running couphng, 

oo 

as{k) ^ as{Q)/[l + 2Poas{Q)lnk/Q] = ^ ag(Q) (-2/?o In A^/Qf-i , (19) 

n=l 

into the right-hand side of Eq. ( [l8[ ) and by integrating term by term. This leads to pertur- 
bative coefficients that grow factorially at high orders: 

Rn+i = n! (1/p) (2/3o/p)" . (20) 

The factorial growth implies that the series for Rpcrtio^s) is not convergent but can be 
interpreted as an asymptotic expansion. One should truncate the series at the order n = 
nmax at which the ratio of two successive terms become of order unity and estimate the 
truncation ambiguity SRpert by the size of the last term that is neglected. The evaluation 
of Umax from Eq. (|^) gives 

(i?„+iag+^(Q))/(i?„ag(Q)) ~ 1 ^ n„,, ~ p/ (2/?oa5(Q)) • (21) 

Using this value of rimax and the QCD expression as{Q) ~ 1/{2(3q IuQ/Aqcd), one even- 
tually obtains an ambiguity of the perturbative component that has a power-like behaviour 
with respect to the large scale Q: 

6Rpert ~ (Rn+iaT'iQ))^^^^^^ ~ e""™'"- ~ {Aqcd/qY . (22) 

The recipe as -^ Ois{k) of Eq. ( p!8| ) to estimate the large-order behaviour of the pertur- 
bative series ca,u be justified by summing Feynman gr aph s in the limit of a large number 



iphaJ, 
iiiyii2 



Nf of flavours tl3 [Nj -^ oo). Further justifications 111111113 follow from the structure of 
soft-gluon resummation formulae (cf. Eq. (|^). 

Two main features of the renormalon approach to power corrections are its predictivity 
and its non-universality. The power p in Eq. (^) is exactly equal to that on the right-hand 
side of Eq. (|lq). Therefore, the approach unambiguously predicts the type of power correc- 
tion by relating it to the dominant low-momentun behaviour of the lowest-order Feynman 
diagrams for any given observable. Nonetheless, the actual size of the power correction is not 
computable unambiguously: the 'true' coefficient in front of the factor 1/Q^ is not related 
in a straightforwar d ar id universal way to that obtained by evaluating the sole lowest-order 
Feynman diagrams 1113. 

5.2 Hadronic event shapes in e^e~ and DIS 

The renormalon predictions on the power behaviour of non-perturbative corrections are in 
agreement IHj with those from OPEs when the latter apply, but can be extended also to 
other quantities. 

Among these quantities, hadronic event shapes in e^e~ annihilation are particularly 
relevant because their measured mean values receive significant non-perturbative contribu- 
tions of the form 1/Q. As shown in Fig. ^, parton level predictions of Monte Carlo event 
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Figure 10: Energy dependence of the mean value of e'^e event shapes. The upper (lower) curves are the 
predictions of Monte Carlo event generators at hadron (parton) level. 



generators fail to describe the data, but a good agreement is obtained by taking into ac- 
count 1/Q effects produced by tlie hadronization models that are built into the Monte C arlo 
programs. The phenomenological success of these hadronization models is well known 1113, 
but their relation with QCD dynamics is still poorly understood on a theoretical basis 

An important result of the renormalon approach is that it theoretically predicts £1^ the 
existence of 1/Q non-perturbative corrections to the mean values of event shapes. As for 
the actual size of these corrections, an attempt to ove rcome non-universality in renormalon 



(DW). The factorial growth of 
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calculations was proposed by Dokshitzer and Webber 
the perturbative coefficients in Eq. (20) is produced by the A;- integration (see Eq. (18)) of 
the perturbative coupling as{k) down to the Landau pole at k = Aqcd- The DW model 
assumes that a meaningful non-perturbative definition of as{k) can be introduced for all 
values of k. Thus the integral 



/ — asik) = ao(w) 
Jo W 



(23) 



exists for all /x/ > and, using Eq. (|^) with p = I, the mean value {R{Q)) of any event 
shape can be written as follows 

(RiQ)) = {RpcrtiasiQ)))NLO + {Rnon-peAQ)) , (24) 

(i?non-pert(Q)) = Oi? Ooifll) (fJ-l/Q) - ttR [ao{fIl)]j^^Q [P-l/Q) + ■■■ ■ (25) 



The first term on the ri gh t-hand side of Eq. ( |24D is the customary perturbative contribution 
evaluated up to NLOB3O, while the power correction (-Rnon-pert(Q)) is expressed in terms of 
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Figure 11: QCD fits for tfie mean value of e^e event shapes. The solid lines correspond to NLO calculations 
combined with power corrections. 



the effective non-perturbative parameter ao(^/)- To avoid double-counting of perturbative 
contributions, one has to consider the NLO expansion [ao(/^/)]AfLo °^ ^^^ integral ( p3D 
and then subtract the second term on the right-hand side of Eq. (^). The coefficient a^ 
depends on the event shape R and is obtained by a direct calculation of the corresponding 
LO Feynman diagrams. 

The DW model does not predict the absolute value of (-Rnon-pcrt(Q)) for each event 
shape R but parametrizes all these power corrections in terms of the single 'universal' 
parameter ao(/^7) and of calculable process-dependent coefficients a/j. The dots on the right- 
hand side of Eq. ( p5[ ) stand for contributions that are more power-suppressed (i.e. {(ii/QY 
with p > 1) and for non-universality corrections to the coefficient a ^ of the dominant power. 



These corrections come from higher-order Feynman diagrams 
integrals of the type 

JO A*/ 



LIS 



and, in particular, involve 



(26) 



If the non-perturbative parameter ao(M/) turns out to be relatively small, the corrections 
in Eq. ([2q ) can be neglected in a first approximation. In phenomenological applications, 
the infrared matching scale /i/ of Eq. ( |2^ ) has to be chosen in the range Q ^ jjli ^ ^qcd- 
Q ^ ^i, because power corrections of higher order have to be negligible and /x/ » Aqcd, 
because as{m) still has to be in the perturbative region. 



The DELPHI Collaboration performed a detailed studyliiS'liiJ of the energy dependence 
of the mean value of e"'"e~ event shapes measured in the centre-of-mass energy range Q = 14- 
172 GeV. The data are well described by NLO perturbative calculationsBjcJ supplemented 
with power corrections (Fig. \\A\j . Using the predictions of the DW model for the thrust T 
and the heavy jet mass M^ (Fig. p|b), a combined fit of as{Mz) and the non-perturbative 
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Figure 12: Distribution of the thrust (Tz) of the current fragmentation region in DIS at different values of 
Q. The histograms are the corresponding NLO predictions. The spectra are multiplied by factors 10", n = 
0, 1, 2, 3 (downwards). 



parameter ao(w) giveslll3 asiMz) = 0.116 ± 0.002(exp.) ± 0.006(th.) and 



aoim = 2 GeV) 



0.534 ±0.012 
0.435 ±0.015 



(from T) , 
(from Mh) 



(27) 



This determination of ^^(M^) is consistent with those (cf. Sect. |3.1| ) obtained from analyses 
of shape distributions in which non-perturbative effects are estimated by using Monte Carlo 
event generators. The value of oq in Eq. (27) suggests that corrections in the DW model can 
be kept under control. Similar results have been obtained by the JADEe3 and ALEPHtH3 
Collaborations. 

Note that the fit of Fig. p| b for the jet mass difference M^ corresponds to an empirical 
parametrization of the power correction simply proportional to 1/Q rather than to the 
predictions of the DW model. In fact, in this case the model gives a poor quantitative 
description of the data. This fact signals the presence of non-universality corrections, which 
are not included in the naive version of the model. 

The data in Fig , lie are obtained by averaging out the event shapes over a restricted 
kinematic range tla that excludes the two-jet region. Note that, in the cases of T and 



Mh, the energy dependence of the data in Fig. |ll| c is stronger than that in Fig. |ll|b. The 
empirical fits of Fig. pTJ c are consistent with a dominant power-behaviour of the type 1/J^ 
This behaviour agrees with that predicted by calculations of renormalon contributions 
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Lepton-nucleon DIS processes are an ideal place to study the Q-dependence of hadronic 
observables because one can vary Q over a wide range in a single experiment, thus reduc- 
ing systematic uncertainties. Event shapes in DIS are analogous to those in e~^e~ an- 
nihilation, but they are usually defined in the current fragmentation region of the Breit 



frame. A study liiitl of DIS event shapes in the momentum region Q= 7-100 GeV has re- 
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Figure 13: Q-dependence of the mean value of DIS event shapes. The sohd hues are QCD fits that include 
NLO predictions and power corrections according to the DW model. 



cently been performed by the HI Collaboration at HERA. As shown in Fig. ^, shape 
distributions becon ie narrow er and jet-like when Q increases. As expected, NLO pertur- 
bative calculations 



i 



123,123,12^ 



describe the data only at sufficiently large values of Q. At 
smaller Q, non-perturbative contributions are sizeable. The Q-dependence of the mean 

together with a combined fit that uses 



value of the sha pp yanabl Rs is shown in Fig. 13 



NLO predictions 
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123 



125 



The fit describes 
(th.) and ao{ni = 2 GeV) = 



'12a and the DW model for power corrections 
the data weU and givesllil as{Mz) = 0.118 ± 0.001(exp.)+[};[]J]^ 

0.491 ± 0.003(exp.)lQ;Q42(th.). These values are consistent with those obtained from e'^e~ 
event shapes. 

Determinations of as from QCD predictions that include power corrections are sum- 
marized at the end of Table 1. Besides the results from the mean values of shape variables, 
I have also reported a preliminary determination IIHj from the longitudinal cross section in 
e^e~ annihilation. This determination is based on the NLO calculation in Ref. liHZI and on 
estimates llHa of power-suppressed terms. 

At present, it is difficult to quantify the theoretical accuracy of these analyses. Uni- 
versality of power corrections is certainly violated in its naive formllHHI. In the DW model, 
the effective non-perturbative parameter is ao ~ 0.5 and, on a phenomenological basis, one 
may conclude that non- universality pffpct s ar p t ypic ally smaller than 50%. Future investi- 



gations on differential event shape spectra 



119 



12C 



13C 



can produce more definite quantitative 
results and give additional insights into the connection between hadronization and power 
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corrections. 



6 Summary of as 



A summary of as determinations is presented in Table [l] I shall limit myself to few 
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comments on the comparison with last year summaries ■ 

The most relevant new result regards a redetermination 11^ of as from i^-nucleon DIS. 
The preliminary result, as{Mz) = 0.119 it 0.005, presented by the CCFR Collaboration at 
the Warsaw Conference (ICHEP96), has indeed been confirmed. This result is based on a re- 
analysis of the CCFR data, due to a new energy calibration of the detector, and supersedes 
a previous and lower determination {as{Mz) = 0.111 it 0.006). The main outcome of this 
is an increased average value of as from DIS measurements. Among DIS determinations, 
only /i-nucleon data from BCDMS still prefer E^ a value of as{Mz) that is (slightly) lower 
than those from e^e~ annihilation. 



Another new result IH3 included in Table |l] is from scaling violation of polarized DIS 
structure functions. In the case of polarized DIS processes, a more accurate determinationEJ 
of as is in principle feasible from the Bjorken sum rule. At present, however, this deter- 
mination suffers from uncertainties t^ due to the extrapolation of the data in the small-x 
region, where polarized structure functions have not yet been measured. 

The entry from pp — > jet -|- X does not refer to an actual measurement, but it illus- 



tratesli^ the potential of jet data from hadron colliders in the determination of as- Detailed 
analyses at j-. he Tevatron experiments are in progress. ^ Th e updated entry from r decay is 
from Ref. 
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The theoretical uncertainty on the value&2a from J/^ and T decays has been 
slightly increased according to estimates H^ of additional contributions from colour-octet 
operators. 

All the other new or updated entries in Table |l| have been discussed in the previous 
sections. In particular, the error on as from global fits to electroweak observables (cf. 
Sect. p.l| ) includes the theoretical uncertainty due to the SM assumptions. 

The values of as, as a function of the energy scale Q at which they are measured, are 



compared with the QCD prediction of a running coupling in Fig. |lj. The energy dependence 
of as is distinct {as{Q) varies by a factor of 3) and in very good agreement with the QCD 
running over two orders of magnitude in Q. It is thus meaningful to evolve all the results 
to as{Mz) according to perturbative QCD. A significant subset of these values is shown in 
Fig. ^. The subset is chosen by considering the most relevant determinations of as from 
each type of process and/or energy range. Using these results, my preferred world average 
determination is 

as{Mz) = 0.119 ± 0.005 . (28) 

The central value corresponds to the weighted average of the results in Fig. |l^ (the average 
does not include the result from fcift mass splitting because of the difficulty in estimating 
uncertainty in lattice calculations cl3). The errors on the most precise determinations of as 
are mainly of a theoretical nature and, hence, not Gaussian and highly correlated among 
themselves. Because of this reason, I have not considered any further reduction of the error 
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Figure 14: Summary of as{Q)- The lines indicate the QCD predictions in NNLO for three different values 
of as (Mz). 



coming from the weighting procedure. The uncertainty quoted in Eq. (^) is equal to the 



smahest of the errors in Fig. 15. 



Different methods, used in recent summaries 
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m 



results similar to that 
are considered in Refs.E^ 



!)• 



to quote the world average of as, give 
Less conservative error estimates {Aas{Mz) = 0.003) 



7 Concluding remarks 

Owing to high-precision experiments and accurate theoretical calculations, perturbative 
QCD is nowadays very well tested in high-energy hadronic processes. An oustanding ex- 
ample of this is the accuracy on the determination of the strong coupling as{Mz) and 
its running. Taken literally, this accuracy implies that in many processes we can control 
strong-interactions dynamics at short distances to a precision better than 5%. The achieved 
reliability of perturbative QCD is extremely valuable to estimate SM backgrounds for new- 
physics signals, although in some cases the insufficient knowledge of the non-perturbative 
parton densities remains a source of sizeable uncertainty. 

Small-x physics deals with a kinematic regime near the borderline between hard and soft 
collisions and can lead to a QCD understanding of the high-energy behaviour of soft hadronic 
interactions. In the very near future substantial progress is expected as a consequence of 
the increasing amount of data and of recent theoretical developments in this field. 

The interplay between perturbative and non-perturbative phenomena is one of the main 
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Figure 15: Summary of measurements of as{Mz)- The band corresponds to the world average as{Mz) = 
0.119 ±0.005. 

open issues in QCD. In the past few years, methods inspired by perturbation theory have 
been developed to control non-perturbative power corrections to perturbative predictions. 
These methods are at present in a stage that is similar to that of LO perturbative calcula- 
tions at the end of the 70's. Forthcoming phenomenological and theoretical studies along 
these lines can play an important role in understanding the hadronization mechanism. 
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Table 1: A summary of measurements of as- 



Process 
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cxp. 
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Order of 
perturb. 


GLS sr 


1.73 
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0.115 ±0.006 


0.005 
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NNLO 


Rr 
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0.002 


0.005 


NNLO 


DIS [polar.] 
DIS [ HERA Fa] 
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Q oi + 0.08 
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0.23 ±0.04 


120 + "•''^'' 

U-J-ZU _ 0.008 

0.120 ±0.010 


+ 0.004 
- 0.005 

0.005 


+ 0.009 
- 0.006 

0.009 


NLO 
NLO 


DIS [ly] 


5.0 


0.215 ±0.016 


0.119 ±0.005 


0.002 


0.004 


NLO 


DIS l/i] 


7.1 


0.180 ±0.014 


0.113 ±0.005 


0.003 


0.004 


NLO 


bb mass sphtting 


4.1 


0.223 ± 0.009 


0.117 ± 0.003 


0.000 


0.003 


LGT 


e+e- [T ± X] 
J/*, T [had. decay] 


4.1 
10.0 


n DOS + 0.045 
U.ZZO _ 0.030 

0.167 ±0.020 


n no + 0010 
u.iiy _ 0.008 

0.113 ±0.010 


0.002 
0.001 


+ 0.010 

- 0.008 

0.010 


NLO 
NLO 


e+e" [ahad] 
e+e" [ev. shapes] 
e+e" [ev. shapes] 
e+e" [ahad] 
e e" [ev. shapes] 
e+e" [ev. shapes] 
e+e~ [ev. shapes] 


10.5 
22 
29 
34.0 
35.0 
44.0 
58.0 


0.20 ±0.06 
0.161 t HI? 
0.160 ±0.012 

0.146 t °:«^J 
0.143 1 °:S^ 

r, ,07 + 0.010 
U.i^^ _ 0.007 

0.132 ±0.008 


0.13 ±0.03 

r, 194 + 0.009 
U.i^4 _ 0.006 

0.131 ±0.010 

U.i^4 _ 0.019 
n 199 + 0.008 
U.iZZ _ 0.006 
r, ,99 + 0.008 
U.l^^ _ 0.006 

0.123 ±0.007 


0.02 
0.005 

0.006 

+ 0.021 

- 0.019 

0.002 
0.003 
0.003 


0.02 

+ 0.008 

- O.ooa 
0.008 

+ 0.008 

+ §:88? 

- 0.005 

0.007 


NNLO 
resum. 
resum. 
NLO 
resum. 
resum. 
resum. 


pp,pp ^ y + X 
pp^ bb + X 
pp -^ jot ± X 


4 
20.0 
120 


0.206 + nS 
0145 + o'oi8 

U.140 _ 0.019 

0.116 ±0.009 


112 + "-"^^ 

^•-'--'-^ - 0.008 

0.113 ±0.011 
0.121 ±0.010 


0.006 

+ 0.007 
- 0.006 

0.008 


+ 0.010 

- 0.005 
+ 0.008 

- 0.009 

0.005 


NLO 
NLO 
NLO 


e+e" [seal, viol.] 


36 


0.147 ±0.014 


0.125 ±0.010 


0.006 


0.008 


NLO 


Z° [e.w. obs.] 
Z° [ev. shapes] 


91.2 
91.2 


120 + "-""^ 

U.lZU _ 0.004 

0.122 ±0.006 


r, 190 + 0.005 
U-J^^U _ 0.004 

0.122 ±0.006 


0.003 
0.001 


+ 0.004 
- 0.002 

0.006 


NNLO 
resum. 


e+e" [ev. shapes] 


133 


0.111 ±0.007 


0.117 ±0.008 


0.004 


0.007 


resum. 


e+e" [ev. shapes] 


161 


0.106 ±0.007 


0.115 ±0.008 


0.004 


0.007 


resum. 


e+e" [ev. shapes] 


172 


0.103 ±0.007 


0.112 ±0.008 


0.004 


0.007 


resum. 


DIS [av. ev. shapes] 
e e"[av. ev. shapes] 
e+e" [aL] 


7- 100 

14 - 172 

91.2 


0.110 ±0.016 


118+ 0.007 
U.iiS _ 0.006 

n lift + o""^ 

U.ilD _ 0.005 

0.110 ±0.016 


0.001 
0.001 
0.013 


+ 0.007 

+ 0.006 
- 0.005 

0.010 


NLO ± p.c. 
NLO ± p.c. 
NLO ± p.c. 
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